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ABSTRACT 

Motivated by unexplained observations of low sulphur abundances in planetary nebulae (PNe) 
and the PG1159 class of post asymptotic giant branch (AGB) stars, we investigate the pos- 
sibility that sulphur may be destroyed by nucleosynthetic processes in low-to-intermediate 
mass stars during stellar evolution. We use a 3 M©, Z = 0.01 evolutionary sequence to exam- 
ine the consequences of high and low reaction rate estimates of neutron captures onto sulphur 
and neighbouring elements. In addition we have also tested high and low rates for the neutron 
producing reactions ^^C(a,n)'''0 and ^^Ne(a,n)^^Mg. We vary the mass width of a partially 
mixed zone (PMZ), which is responsible for the formation of a '^C pocket and is the site of 
the '^C(a,n)'^0 neutron source. We test PMZ masses from zero up to an extreme upper limit 
of the entire He-intershell mass at 10^^ M©. We find that the alternative reaction rates and 
variations to the partially mixed zone have almost no effect on surface sulphur abundances and 
do not reproduce the anomaly. To understand the effect of initial mass on our conclusions, 1.8 
M0 and 6 M^, evolutionary sequences are also tested with similar results for sulphur abun- 
dances. We are able to set a constraint on the size of the PMZ, as PMZ sizes that are greater 
than half of the He-intershell mass (in the 3 Mp, model) are excluded by comparison with neon 
abundances in planetary nebulae. We compare the 1.8 M© model's intershell abundances with 
observations of PG1159-035, whose surface abundances are thought to reflect the intershell 
composition of a progenitor AGB star. We find general agreement between the patterns of 
F, Ne, Si, P, and Fe abundances and a very large discrepancy for sulphur where our model 
predicts abundances that are 30-40 times higher than is observed in the star. 

Key words: ISM: abundances; planetary nebulae: general; stars: evolution, AGB and post- 
AGB. 
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1 INTRODUCTION 

After leaving the asymptotic giant branch (AGB), a post-AGB star 
may evolve to high temperatures (> 30,000 K) on the timescale 
required to ionise the surrounding shell of ejected material and be- 
come visible as a planetary nebula (PN). As PNe are composed of 
envelope material from a progenitor AGB star, measurements of PN 
chemical abundances provide a way to test the predictions of AGB 
nucleosynthesis models (e.g., 
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[Marigo et al.|[2()03 

[2009i|Pottasch & Bemard-Salas|2010[|Karakas & Lt igaro 20101. 

Many planetary nebulae (PNe) with approximately solar oxy- 
gen abundance (±0.4 dex) have been found to have sulphur deple- 
tions of between 0.1 and 0.6 dex relative to the Sun ( [Marigo et al.| 
|2003| l. A more detailed investigation with a larger sample of 85 
PNe by |Henry et al.|pb04^ discovered that sulphur abundances in 
PNe are systematically lower than HII regions at the same metal- 
licity, where metallicity in PNe is measured indirectly through the 
oxygen abundance. Specifically, [Henry et al.| ( [2004[ l showed that 
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the abundance trends between PNe and HII regions are co-linear 
in the Ne-0, Cl-O, and Ar-0 planes, but are separated in the S-0 
plane, in which the trend-line of PNe is located below that of HII 
regions by 0.3 dex. The co-linear trends between Ne, CI, Ar, and 
O, but not S single out sulphur as the anomalous element, and this 
has been labelled the 'sulphur anomaly'. The sulphur anomaly has 
been independently confirmed by the observations of |Milingo et al.[ 
pOTOt . 

[Henry et al.[ ( [2004| argued that because AGB models do not 
predict significant depletion of sulphur and sulphur does not read- 
ily condense into dust grains, the most likely cause of the sulphur 
anomaly is a failure to correctly account for sulphur in the highly 
ionised S+"' state through the use of an lonisation Correction Fac- 
tor (ICF) and measurements of S+' and S+^ abundances. However, 
infrared observations of PNe (e.g., [Bemard-Salas et al.|2008[ l have 
directly measured S+"' abundances using the [S IV] emission line at 
10.5 /im. This was done without the need for an ICF and these ob- 
servations show that the sulphur anomaly still exists and is in need 
of explanation. 

If the observed low gas phase abundance of sulphur in PNe 
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relative to the interstellar medium (as sampled by HII regions) re- 
flects a decrease in elemental sulphur between the birth composi- 
tion of a star and its surface layers at the final phases of stellar evo- 
lution, then an attractive solution would be to identify a nucleosyn- 
thetic process that is able to destroy sulphur during the interven- 
ing stages. The progenitors of planetary nebulae are typically low- 
mass stars (e.g., 1.0 to 2.5 MQ; |Pottasch et al.|2011^ , which evolve 
through the AGB phase and experience nucleosynthesis through H 
and He burning and the slow neutron capture process (^-process) 
(e.g., |Gallino et al.|1998[|Busso et al.|i999) . During the AGB, nu- 
cleosynthesis products are periodically dredged up into the convec- 
tive hydrogen-rich envelope as part of the thermal pulse cycle, so 
a depletion of sulphur in the He-intershell of an AGB star would 
result in a (smaller) depletion of sulphur at the stellar surface. A re- 
view of AGB evolution and modelling is given by |Herwig| ( (2005l >. 

[Henry et al.H2012) provides an update on the status of the sul- 
phur anomaly and discusses the still-viable explanations, includ- 
ing gas phase depletion due to dust or molecule formation, and the 
nuclear processing in AGB stars. Although it was argued that the 
sulphur anomaly is inconsistent with the predictions of existing nu- 
cleosynthesis models (e.g., |Karakas|2010K there has not been an 
investigation into how modelling uncertainties such as nuclear re- 
action rates and the treatment of mixing affect predictions of sur- 
face sulphur abundances. 

There is a separate physical site with unexplained sulphur de- 
pletion in stars of the type PG1159. PG1159 stars are extremely 
hot (75,000-200,000 K) post- AGB stars that are hydrogen deficient 
and helium rich, likely because of a late or very late helium shell 
flash that has consumed their remaining hydrogen envelope and 
exposed He-intershell material to the stellar surface ( |van Winckel| 
|2003[[Wemer et aL|l 991 1. Although a very late thermal pulse and 
hydrogen ingestion episode may lead to some additional light el- 
ement and i-process nucleosynthesis after the AGB ( [Herwig et al.| 
|201 1[ |Stancliffe et al.|2011| l, the resulting surface abundances are 
expected to largely reflect the intershell composition at the end of 
the AGB phase. With intershell matter at their surfaces, PG1159 
stars provide a test of nucleosynthesis models that is relatively free 
of the uncertainties related to dredge-up efficiency that affect the 
surface abundances of AGB stars. 

[Werner & Herwigl ( [2006| l report that PGl 159 stars have highly 
scattered and generally low sulphur abundances ranging from 0.01- 
1 times solar, while the models of Herwig show He-intershell sul- 
phur abundances at the end of the AGB that are 0.6-0.9 times so- 
lar They suggest that a study is needed to understand how the un- 
certainties of neutron capture cross sections affect intershell abun- 
dances. |^femer£Far|([2009j interprets the discrepancy between low 
sulphur observations of PGl 159 stars and sulphur preserving theo- 
retical models as a failure of stellar modelling. If the current mod- 
els' failure to reproduce the sulphur anomaly is an indication that 
our understanding of stellar nucleosynthesis is in need of refine- 
ment, then a solution to the sulphur problem may lead to a better 
understanding of other aspects of stellar nucleosynthesis, such as 
the mixing near convective boundaries and nuclear reaction rates. 

Sulphur has a significant nuclear charge (Z = 16) that prevents 
it from strongly participating in charged-particle reactions (e.g., p- 
and a-capture) at the temperatures and densities of AGB stellar 
interiors. For sulphur destruction, the much more likely pathway 
is neutron capture on S to produce the unstable isotopes ■'^S and 
^^S, which spontaneously decay via P~ to ^^Cl and ■'^Cl, respec- 
tively. Although neutron-capture reactions with S are included in 
existing models, there is disagreement over the rates of these re- 
actions that becomes particularly significant at temperatures below 
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Figure 1. Reduced reaction rates for '^S(n,y)'^S as a function of tempera- 
ture from several sources in the JINA Reaclib database jCyburt et al.|2010") . 
Source labels are defined in Sectionim 



lO' K. The disparity is evident in Figure[l] which shows the rate of 
the ^^S(n,7)-'^S reaction as a function of temperature from several 
sources in the JINA Reaclib database ( [Cyburt et al.|[2010| . Most 
neutron captures at the top of the He-intershell take place at a tem- 
perature of the order of 10^ K, and at this temperature the indepen- 
dently predicted rates disagree by up to factor of 10. At tempera- 
tures less than lO' K, the neutron capture rates diverge rapidly and 
differ by over three orders of magnitude. 

In this paper, we aim to determine whether the sulphur 
anomaly in PNe is the consequence of the nucleosynthetic pro- 
cesses in PN-progenitor AGB stars. To do this, we calculate models 
that span the range of relevant uncertainties in modelling low-mass 
stars - neutron-capture reaction rates, neutron-producing reaction 
rates, and partial mixing zone profiles (which determine the size 
of the '■'C-pockets) and compare their surface abundances to PNe 
observational data. We will examine the intershell abundances of 
our models in comparison with PGl 159 observations and consider 
the significance that extra mixing (e.g., convective overshoot and 
rotation) could have reproducing the sulphur anomaly in PNe. 



2 NUMERICAL METHOD & MODELS 

We evolve our stellar evolutionary sequences from the zero-age 
main sequence to the tip of the AGB with the Mount Stromlo Stel- 
lar Structure Program, which has been updated to include C- and N- 
rich low temperature opacity tables from [Lederer & Aringer[ ( |2009| l, 
as described in [Karakas et al.[ ( |2010| l and references therein. We use 
Reimers' formula jReimers|I975 > with the parameter Tj = 0.4 for 
mass loss during the first red giant branch and the [Vassiliadis &[ 
[Wood[(T993| prescription on the AGB. 

The evolutionary code operates on a minimal set of nuclides 
that are involved in reactions that are highly exothermic (the pp- 
chains, CNO cycle, triple-O!, and '^C(a,y)'*0 reactions) and hence 
affect the stellar structure. The structure model generated by the 
evolutionary code is used as input to a post-process nucleosynthe- 
sis code. With time- and mass-dependent variables such as temper- 
ature, density, and the locations of convective boundaries defined 
in the structure model, the nucleosynthesis code recalculates abun- 
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dances for a detailed network with time-dependent diffusive mixing 
for all convective zones jCannon|1993[ l. 

As the energy generation in the He-flash convective zone 
is completely dominated by the triple-a reaction included in the 
structure model, we are able to modify the rates of weakly ener- 
getic or endothermic a- and n-capture reactions and add a partially 
mixed zone (Section [2T2| in the nucleosynthesis post-process with- 
out having to recalculate the stellar structure. In this work, our de- 
tailed nuclear network consists of 125 species, which include many 
isotopes of P, S, and CI to precisely account for the neutron-capture 
and j3 -decay reactions around sulphur. 

Our chosen initial mass of 3.0 Mq is near the upper end of PN 
progenitor masses. As shown in the results of |Karakas & Lattanzio| 
( |2007^ , 3.0 Mq models experience a greater number of thermal 
pulses and third dredge-up than lower mass models, so this choice 
will exaggerate the effect of any possible sulphur depletion on sur- 
face abundances. 

For our initial composition, we scale the solar abundances of 
[Asplund et"aL] ( |2009| (Zg = 0.0142) such that the models' metal- 
licity is Z = O.7O4Z0 = 0.01. The models' final metallicity will 
be roughly a factor of two larger than the initial value (mostly due 
to the dredge-up of primary '"C during the thermally pulsing AGB 
phase), so the model will lie just above the centre of the PNe metal- 
licity range of 0.3 Zg to 2Z0 ^Sterling & Dinerstein|2008] l. 



2.1 Changes to Reaction Rates 

To explore the effect of rate uncertainties for neutron-capture re- 
actions with sulphur and its nuclear neighbours, we select the two 
sources in the JINA Reaclib database that predict the highest and 
lowest rates for ^^S(n,7)^^S around the intershell temper ature of 
10^ K. Our standard rate case is the Reaclib V0.5 release by Cyburt 



|et al.| ( (2bl0| l, which, for this reaction includes experimental esti 
mates from the KADoNiS database ( [Dillmann et al.|2006| l labelled 
'ka02'. The source that predicts the lowest rate, 'kd02' is very sim- 
ilar to 'ka02', except that fitting formulae have been adjusted to 
maintain accuracy at low temperatures. The highest rate source, 
'ths8' is comprised of theoretical estimates by Thomas Rauscher 
that were included as part of the REACLIB VI.O release jCytart] 
|et al.|2010| l. The cuiTent ReaclibV2.0 release has adopted 'kd02' 
rates for the ^^S(n,y)-'"'S reaction. Also included in Figure[I]are the 
experimental rates of Bao & Kappeler ( 1987 1 (labelled 'baka'), and 
the statistical model calculations by Rauch & Thielemann| f2000^ 
('rath'), however, these rates have not been adopted in this study. 



2.2 Partial Mixing Zone 

The free neutrons available for the .s-process in low-mass stars 
are primarily produced by '^ C burning under radiative conditions 
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via th e '^C(a,n)'^0 reaction i Straniero et al.|l995 
19981. Producing the seed '^C nuclei in stellar models requires 



the existence of a layer at the top of the '^C-rich intershell in 
which protons are 'partially mixed' down from the envelope, thus 
enabling the CN cycle reaction i2c(p,y)i4N(j3 + )'3C. The rmxmg 
process cannot be too efficient, or else the newly-created ^^C nuclei 
will be destroyed by further proton capture to make ^^N, which is 
a neutron poison, i.e., its large neutron-capture cross section sig- 
nificantly reduces the number of free neutrons available for the s- 
process. The physical mechanism behind the formation of a par- 
tially mixed zone (PMZ) is still a mystery, although some of the 
more likely possibilities include convective overshooting jHerwig| 
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Figure 2. Upper panel: The 3 M j model proton profile immediately after 
inserting a 10^' M,;; partial mixing zone. Lower panel: As the envelope 
convection zone (shaded) retreats outwards in mass, proton capture reac- 
tions result in adjacent pockets of '^^C and '''N forming near the top of the 
intershell. Y denotes the molar fraction, equal to (mass fraction) / (atomic 
mass). 
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|Cristallo et al.|2004) l, rotational mixing fHerwig & Langer| 



or gravity-wave driven mixing (Denissenkov & Tout 2003 
Some models in the literature insert a '^C pocket directly at 
each thermal pulse, using a profile such as the |Gallino et al.|(l998| ) 
standard (ST) case, which has a '^C pocket mass of 5 X 10"^ Mq . 
[Herwig & Langerj ( (200T] l use a diffusive convective overshoot at 
the bottom of the envelope convection boundary with the parameter 
/ = 0.016 and find a '^C pocket width (where the '^C mass fraction 
is above lO"'*) of about 2 x 10"^ M© in a 3 Mg model. 

To experiment with different '^C pocket masses, the studies 
o flArlandini et al.|(fT999l l and |Bisterzo et al.J pOTOl [20T2) insert 
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been scaled in C (and ''*N) abundance. We instead scale the width 
(in mass coordinate) of an inserted proton profile, which not only 
controls the total mass of protons inserted (and the mass of the 
resulting '^C and '^N pockets) but also changes the radial position 
and extent over which the resulting neutron-captures take place. For 
a comparison involving both of these treatments of the '^C pocket, 
see the detailed discussion in |Lugaro et al.|pbl2^ . 

Figure [2] shows the '■'C and ''^N pockets that form as a result 
of an exponential profile of protons inserted below the envelope 



( 1998^ ST profile '^C (and ^^N) pocket that has 
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convection zone. The protons are inserted identically after every 
thermal pulse with third dredge-up, at the time when the envelope 
convection zone reaches its deepest extent during a third dredge- 
up episode. The proton profile matches the envelope abundance at 
top of the PMZ and decreases exponentially to a mass fraction of 
10^^ across a mass interval we refer to as the 'PMZ mass'. This 
method is described in mo re detail in | Lugaro et aL| (| 2004^ and is 
very similar to that used by |Goriely & Mowlav i (20001. 

To explore the results of exaggerated neutron capture nu- 
cleosynthesis, we test 3 Mp, models with PMZ masses of 
(1,5, and 10) x lO^-' M^ in addition to a model with no partial 
mixing zone. Unless otherwise stated, we use our standard PMZ 
with a mass of 1 x 10"^ Mg. The PMZ mass of 10 x 10"^ Mg is 
included as an extreme upper limit for for the 3 Mg model, as this 
profile spans the entire He-intershell and the partial mixing is not 
expected to penetrate into the degenerate C-O core. 



3 MODEL RESULTS 

3.1 New Stellar Evolutionary Models 

The 1.8 Mg model has previously been described in |Karakas et al.| 
( |2010 i We present new 3 Mg and 6 Mg sequences at Z = 0.01 
computed with Mount Stromlo evolutionary code. 

Table[T]and Table|2]list the structural parameters of each ther- 
mal pulse (TP) in the 3 Mg and 6 Mg models, respectively. The 
parameters include the total mass (Mtot) and the H-exhausted core 
mass (Mh) measured at the beginning of the TP, the mass dredged- 
up into the envelope (AMqup) and the dredge-up efficiency pa- 
rameter (A = AMoup/'^^H' where A denotes the change between 
the previous and current TP) measured after the TP. Also included 
are the maximum temperatures of the He- and H-burning shells 
(THe-shelb TH-shell) during the TP. In the interpulse period (Tip), 
we sample the maximum temperature at the base of the convec- 
tive envelope (Tbce)> the minimum luminosity (Lmax), the max- 
imum stellar radius (Riiiax)> the minimum bolometric magnitude 
(Mboi), the minimum effective temperature (Tgff), and the carbon- 
to-oxygen number ratio (C/0). 

In our 3 Mg model, the temperature at the base of the con- 
vective envelope never exceeds 5.6 x 10^ K, so there is no pro- 
ton capture nucleosynthesis in the envelope (hot-bottom burning), 
which requires temperatures above about (40-50) x 10^ K. With- 
out active CN-cycling in the envelope, third dredge-up raises the 
surface C/0 ratio, eventually causing the model to become carbon- 
rich (C/0 >1) after the Uth thermal pulse. 

Envelope temperatures in the 6 M t, model are in the range 
(31-85) X 10^ K, high enough that hot-bottom burning is active 
and prevents the surface from becoming carbon-rich up to the last 
(40th) thermal pulse, at which time the C/O ratio is 0.95. At the end 
of the AGB, the C/0 ratio has increased to 1.16. 

The maximum temperature in the helium burning shell gen- 
erally increases with thermal pulse number in the 3 Mg model, 
reaching a temperature of about 300 x 10^ K after the 15th ther- 
mal pulse. Above about 300 x 10^ K, the ^^Ne(a,n)^^Mg reaction 
starts to become active, producing a burst of neutrons at the base of 
the intershell during convective thermal pulses. The ^^Ne neutron 
source operates in addition to the '■'C source, so this indicates that 
the intershell matter will be subject to additional neutron captures 
over the last few thermal pulses. 

The 3 Mg model experiences a total of 22 thermal pulses be- 
fore the AGB phase is terminated by mass loss at a final core mass 
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Figure 3. Mass of the He-exhausted (dashed) and H-exhausted (solid) cores 
during the thermally pulsing AGB phase as a function of time in the 3 
model. 

of 0.68 Mg (Figure [3j. At this time, the total amount of intershell 
matter that has been dredged up into the envelope is 0.120 Mg. 
With an envelope mass of about 2 Mg for most of the AGB, there 
is enough material dredged up that the results of neutron capture 
nucleosynthesis in the He-intershell will be evident in the surface 
abundances of the model. This amount of dredge-up is a factor of 
a few compared to our evolutionary models with initial masses of 
1.8 Mg and 6 Mg, which dredge up a total of 0.041 Mg and 0.082 
M;g of intershell matter, respectively. 

The cumulative dredge-up quantity of 0.12 Mp^ in our 3 
Mg, Z = 0.01 model is comparable to other models of the same 
mass and similar metallicity. [Lugaro et al.| (2003) compare 3 Mg, 
Z = 0.02 models (the same mass but higher metallicity than our 
model) that have been computed independently with the Mount 
Stromlo stellar structure program, FRANEC ( |Gallino et al.|1998| l, 
and EVOL (as used by |Herwig|2000[ with hydrodynamic overshoot 
included at all convective boundaries), finding cumulative dredge- 
up quantities of 0.08 Mg, 0.044 Mg, and 0.10 Mg, respectively. 
A 3 Mg, Z = 0.02 model calculated with the Cambridge STARS 
code, which computes mixing and burning as a single step, dredges 
up a total of 0.13 Mp^ after 20 thermal pulses in [Stancliffe et al.] 
p004t . 

3.2 Nucleosynthesis Model Results 

Table [3] presents the final surface abundance results of the 3 Mg 
model with standard ('ka02') and alternative neutron capture rates 
('kd02' and 'ths8') and several partial mixing zone sizes (0, 1, 5, 
and 10 X lO^-' Mg), as illustrated in FigureUJ The alternative neu- 
tron capture rates are seen to have little effect (< 0.001 dex) on 
the surface abundance of sulphur at the end of the AGB phase. The 
variations to the partial mixing zone size also leave sulphur abun- 
dances virtually unchanged. In each of our cases, the surface abun- 
dance of sulphur is not depleted but instead shows a very slight 
increase (< 0.03 dex) from pre-main sequence to the end of the 
AGB (FigurepV 

In Figurep] we show surface abundance evolution in the S vs. 
O plane for 3 Mg models with a range of PMZ sizes, as well as 
the 1.8 and 6 Mg models for comparison with observational PNe 
abundances from |Pottasch & Bemard-SalasH2010) . Although it is 
unlikely for a 6 Mg star to produce a detectable planetary nebula, 
the 6 Mg abundances are included to show the extent to which our 
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Table 1. Stnictural and dredge-up parameters of the 3 M,.,, Z = 0.01 sequence for each thermal pulse cycle during the AGB. Columns are 
described in Section lTTI 



TP Miot Mh AMdUP ^ Tae-shcll Tn-shell Tbcc Tip Lmax Rmax Mbol Tcff C/O 

# [Mr,,] [Mr^,] [10^3 M,,] [10* K] [10* K] [10* K] [10^ yr] [10^ Lq] [Rq] [mag] [K] 



1 


2.99 


0.607 


0.00 


0.00 


181.6 


53.0 


2.3 


0.00 


3.5 


146 


-4.15 


3797 


0.36 


2 


2.99 


0.610 


0.00 


0.00 


214.1 


57.3 


2.7 


58.27 


5.0 


188 


-4.52 


3775 


0.36 


3 


2.99 


0.614 


0.00 


0.00 


226.3 


58.9 


2.8 


77.52 


5.6 


204 


-4.65 


3749 


0.37 


4 


2.99 


0.619 


0.00 


0.00 


238.2 


60.7 


3.0 


84.12 


6.3 


223 


-4.78 


3715 


0.37 


5 


2.99 


0.625 


0.00 


0.00 


245.6 


62.0 


3.2 


85.53 


6.9 


238 


-4.87 


3682 


0.37 


6 


2.99 


0.631 


0.90 


0.14 


253.0 


63.1 


3.3 


83.52 


7.4 


251 


-4.95 


3580 


0.37 
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8 


2.99 


0.642 


3.71 


0.47 


269.7 


65.4 


3.8 


81.83 


8.5 


278 


-5.10 


3474 


0.48 


9 


2.99 


0.647 


5.24 


0.61 


273.0 


66.2 


4.0 


81.51 


9.0 


290 


-5.16 


3437 


0.63 


10 


2.99 


0.651 


6.57 


0.70 


278.9 


66.8 


4.3 


82.36 


9.5 


303 


-5.22 


3410 


0.84 


11 


2.99 


0.655 


7.32 


0.72 


283.1 


67.3 


4.6 


82.72 


9.9 


313 


-5.27 


3387 


1.07 


12 


2.99 


0.658 


8.01 


0.76 


287.2 


67.6 


5.0 


81.37 


10.3 


318 


-5.31 


3377 


1.30 


13 


2.98 


0.661 


8.49 


0.77 


290.8 


67.8 


4.9 


80.67 


10.7 


352 


-5.35 


3258 


1.53 


14 


2.98 


0.664 


8.77 


0.78 


294.0 


68.0 


5.1 


78.92 


11.0 


371 


-5.38 


3198 


1.78 


15 


2.98 


0.666 


8.81 


0.78 


296.5 


68.2 


5.2 


76.91 


11.3 


391 


-5.41 


3145 


2.03 


16 


2.98 


0.669 


9.00 


0.80 


298.6 


68.3 


5.3 


74.18 


11.6 


409 


-5.44 


3097 


2.27 


17 


2.97 


0.671 


8.89 


0.78 


300.8 


68.5 


0.0 


72.56 


11.8 


426 


-5.46 


3055 


2.50 


18 


2.96 


0.673 


8.83 


0.80 


301.7 


68.6 


5.5 


69.27 


12.1 


444 


-5.48 


3019 


2.74 


19 


2.92 


0.676 


8.99 


0.81 


302.8 


68.7 


5.6 


67.58 


12.3 


460 


-5.50 


2983 


2.98 


20 


2.83 


0.678 


8.98 


0.81 


304.9 


68.8 


5.6 


65.74 


12.5 


482 


-5.52 


2948 


3.23 


21 


2.17 


0.680 


7.85 


0.71 


305.8 


68.8 


5.4 


63.95 


12.6 


581 


-5.53 


2901 


3.48 


22 


1.23 


0.682 


7.09 


0.72 


303.6 


68.4 


3.4 


58.91 


12.5 


843 


-5.52 


2633 


3.79 
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Figure 4. Surface sulphur abundance as a function of thermal pulse number 
during the AGB phase in the 3 M ;,, Z = 0.01 model with a PMZ size of 
10"^ Mq and three neutron-capture rate sources. 



results are dependent on tlie stellar initial mass. An enhancement 
in O by up to 0.2 dex is seen in the 1.8 and 3 Mp, models, but this 
is not enough to explain the observational trend of PNe. Figure |6] 
showing the results in the S versus Ar plane confirms that neither S 
nor Ar are significantly processed during low-mass evolution, and 
that none of the models can account for the trend of low S abun- 
dances in PNe. 
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Figure 5. Surface abundance results in the S versus O plane for 3 M j mod- 
els with PMZ sizes of (1, 5, and 10)x 10"^ Mq, the 1.8Mq model with a 
2 X 10^3 PMZ, and the 6 Mq with no PMZ. Included for comparison are 
the PNe observational data of Pottasch & Bernard-Salas'('2010 l and the in- 
terstellar medium (ISM) trend of Milingo et al. 1 2010 1 from observations of 
HIT regions and blue compact galaxies. 



The sulphur abundance not only changes very little at the sur- 
face, but remains relatively constant throughout the entire stellar 
interior at the end of the AGB. The abundance discontinuity in 
(Figure [7| at the intershell-envelope boundary shows the results 
of a small ■'^S depletion in the intershell by conversion into more 
neutron-rich isotopes, indicated by increases in ^■'S, ■''*S, ^^S, and 
■'^S abundances. The increase in ■'^S abundance indicates that neu- 
tron captures onto unstable ^^S nuclei occur on a timescale com- 
parable to its /3+-decay mean lifetime of 126.3 days (Audi et al 
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Table 2. Stnictural and dredge-up parameters of the 6 M,.,, Z = 0.01 sequence for each thermal pulse cycle during the AGB. Columns are 
described in Section lTTI 



TP Miot Mh AMdUP ^ Tae-shcll Tn-shell Tbcc Tip Lmax Rmax Mbol Tcff C/O 

# [Mr,,] [Mr^,] [10^3 M,,] [10* K] [10* K] [10* K] [10^ yr] [10^ Lq] [Rq] [mag] [K] 



1 


5.95 


0.911 


0.00 


0.00 


236.7 


76.9 


31.4 


0.00 


23.1 


436 


-6.19 


3518 


0.37 


2 


5.95 


0.912 


0.07 


0.09 


247.3 


78.1 


38.6 


3.23 


24.0 


449 


-6.23 


3494 


0.37 


3 


5.95 


0.913 


0.28 


0.30 


256.5 


79.3 


46.9 


3.36 


25.0 


462 


-6.27 


3486 


0.37 


4 


5.95 


0.914 


0.54 


0.51 


264.9 


80.4 


56.0 


3.44 


26.2 


478 


-6.32 


3475 


0.37 


5 


5.95 


0.914 


0.71 


0.60 


272.7 


81.4 


62.7 


3.53 


27.6 


498 


-6.38 


3465 


0.38 


6 


5.95 


0.915 


0.93 


0.73 


279.6 


82.2 


66.6 


3.65 


29.0 


519 


-6.44 


3457 


0.38 


7 


5.95 


0.915 


1.10 


0.79 


286.7 


83.0 


69.4 


3.77 


30.4 


537 


-6.49 


3449 


0.38 


8 


5.95 


0.916 


1.26 


0.84 


293.0 


83.7 


71.6 


3.92 


31.7 


554 


-6.53 


3446 


0.35 


9 


5.95 


0.916 


1.41 


0.87 


299.2 


84.4 


73.5 


4.10 


32.7 


569 


-6.57 


3426 


0.31 


10 


5.95 


0.917 


1.56 


0.90 


304.8 


85.0 


75.3 


4.28 


33.6 


580 


-6.60 


3398 


0.25 


11 


5.95 


0.917 


1.56 


0.84 


310.3 


85.7 


77.0 


4.49 


34.3 


590 


-6.62 


3376 


0.18 


12 


5.95 


0.917 


1.65 


0.88 


313.8 


86.2 


78.5 


4.55 


34.7 


596 


-6.63 


3347 


0.13 


13 


5.94 


0.917 


1.83 


0.90 


319.6 


86.9 


80.0 


4.85 


35.4 


603 


-6.65 


3344 


0.09 


14 


5.94 


0.918 


1.95 


0.93 


323.6 


87.5 


81.2 


4.99 


35.7 


609 


-6.66 


3326 


0.07 


15 


5.94 


0.918 


2.06 


0.93 


328.0 


88.0 


82.1 


5.19 


36.3 


616 


-6.68 


3319 


0.06 


16 


5.93 


0.918 


2.11 


0.92 


330.6 


88.3 


82.8 


5.30 


36.9 


625 


-6.70 


3308 


0.06 


17 


5.93 


0.918 


2.19 


0.93 


334.7 


88.6 


83.2 


5.42 


37.4 


631 


-6.71 


3294 


0.06 


18 


5.92 


0.919 


2.25 


0.93 


338.0 


88.9 


83.6 


5.52 


37.9 


638 


-6.73 


3288 


0.06 


19 


5.91 


0.919 


2.30 


0.93 


340.4 


89.1 


83.9 


5.62 


38.4 


644 


-6.74 


3278 


0.06 


20 


5.90 


0.919 


2.36 


0.94 


342.8 


89.3 


84.1 


5.69 


38.8 


651 


-6.75 


3269 


0.06 


21 


5.88 


0.919 


2.41 


0.94 


344.7 


89.4 


84.4 


5.79 


39.2 


657 


-6.76 


3265 


0.06 


22 


5.86 


0.920 


2.21 


0.84 


346.1 


89.5 


84.5 


5.86 


39.6 


663 


-6.77 


3257 


0.06 


23 


5.84 


0.920 


2.22 


0.84 


347.4 


89.6 


84.6 


5.92 


39.9 


668 


-6.78 


3251 


0.06 


24 


5.80 


0.920 


2.54 


0.95 


347.5 


89.7 


84.7 


5.94 


40.2 


673 


-6.79 


3244 


0.07 


Z.J 


J. Ij 


u.yzu 


Z.JO 


u.y4 


340. J 


QO C 

oy.o 


(S4.JS 


O.Uu 


4U.4 


/o 


^ on 
-D.oU 


jZj I 


U.U/ 


26 


5.69 


0.920 


2.57 


0.94 


350.8 


89.8 


84.8 


6.06 


40.5 


682 


-6.80 


3228 


0.07 


27 


5.60 


0.921 


2.35 


0.86 


354.2 


89.8 


84.7 


6.05 


40.5 


686 


-6.80 


3216 


0.07 


28 


5.46 


0.921 


2.45 


0.89 


350.2 


89.7 


84.5 


6.04 


40.4 


690 


-6.80 


3205 


0.08 


29 


5.17 


0.921 


2.63 


0.95 


352.0 


89.6 


84.1 


6.10 


40.1 


695 


-6.79 


3193 


0.08 


30 


4.85 


0.921 


2.67 


0.94 


352.3 


89.4 


83.3 


6.25 


39.3 


700 


-6.77 


3169 


0.08 


31 


4.54 


0.921 


2.69 


0.94 


348.0 


89.0 


82.4 


6.27 


38.0 


704 


-6.73 


3137 


0.09 


32 


4.23 


0.921 


2.72 


0.94 


351.0 


88.7 


81.2 


6.32 


36.7 


709 


-6.69 


3104 


0.11 


33 


3.94 


0.922 


2.75 


0.95 


347.9 


88.2 


79.5 


6.38 


35.4 


712 


-6.65 


3072 


0.13 


34 


3.64 


0.922 


2.80 


0.95 


354.1 


87.8 


77.4 


6.51 


34.0 


715 


-6.61 


3042 


0.15 


35 


3.36 


0.922 


2.86 


0.95 


351.3 


87.3 


74.3 


6.65 


32.6 


717 


-6.56 


3011 


0.20 


36 


3.08 


0.922 


2.90 


0.94 


352.5 


86.7 


69.1 


6.91 


31.1 


720 


-6.51 


2981 


0.29 


37 


2.79 


0.923 


2.97 


0.95 


354.6 


86.0 


57.5 


7.24 


29.6 


730 


-6.46 


2944 


0.41 


38 


2.50 


0.923 


3.03 


0.95 


352.9 


85.3 


39.2 


7.67 


28.5 


754 


-6.42 


2901 


0.56 


39 


2.21 


0.923 


3.05 


0.93 


357.7 


84.7 


22.3 


7.99 


27.7 


783 


-6.39 


2856 


0.74 


40 


1.92 


0.923 


3.11 


0.97 


359.2 


84.1 


6.7 


7.92 


27.2 


809 


-6.37 


2817 


0.95 



Table 3. Log E* surface elemental abundances, carbon-to-oxygen ratio (by number) and '^C/'^C isotopic ratios 
at the end of the AGB in our 3 M, . „ Z = 0.01 model. 





Ne 


Mg 


Si 


S 


P 


CI 


Ar 


C/O 


12C/13C 


scaled solar initial (Z = 0.01) 


7.802 


7.472 


7.382 


6.992 


5.242 


5.062 


6.272 


0.550 


89.4391 


m3 zO 1 -standard-pmz 1 
m3z01-kd02-pmzl 
m3zO 1 -ths8-pmz 1 


8.266 
8.267 
8.265 


7.546 
7.547 
7.545 


7.416 
7.416 
7.417 


7.024 
7.023 
7.022 


5.326 
5.318 
5.313 


5.133 
5.134 
5.152 


6.302 
6.302 
6.302 


4.279 
4.279 
4.280 


305.146 
305.329 
305.278 


m3z01-ths8-pmzO 
m3zO 1 -ths8-pmz5 
m3z01-ths8-pmzlO 


8.167 
8.502 
8.627 


7.531 
7.602 
7.640 


7.417 
7.418 
7.419 


7.023 
7.020 
7.021 


5.280 
5.391 
5.407 


5.136 
5.179 
5.186 


6.302 
6.301 
6.300 


4.549 
3.467 
2.874 


318.118 
259.238 
219.838 



*loge(X) = log(A'x/A'«) + 12 
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Figure 6. Same as Figurejsjbut for S versus Ar. 
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Figure 8. Same 


as Figure[5]but for Ar versus Ne. 
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Figure 7. Abundances in the He-intershell after the second last thermal 
pulse of the 3 M .,, Z = 0.01 model with 'ths8' rates and a PMZ mass of 
1 X 10^3 M,;v The shaded regions indicate convective zones. 



|2003i The sulphur depletion (by 22% or 0.1 dex) in the intershell 
is too small to resolve the sulphur anomaly in PNe, which requires 
sulphur depletions of typically 0.3 dex (and up to 0.6 dex) in the 
hydrogen-rich envelope. 

From the results in Table[3] we see that the neon surface abun- 
dances increase significantly with increases to the mass of the par- 
tial mixing zone. The elemental increase in neon is due to ^^Ne 
production from primary '^N in convective pulses via the reaction 
chain '■*N(a,7)'^F(jS + )'**0(a,r)22Ne. The size of the PMZ corre- 
lates with the size of the resulting '^N pocket, which adds to the 
quantity of '^N available for the production of ^^Ne, as shown in 
Figure [5] 

The position of our models with large PMZ sizes in the Ar ver- 
sus Ne plane is far from the observational PNe data (Figure [8j and 
from this we conclude that the partial mixing zone widths larger 
than 5 x 10"^ Mg (50% of the He-intershell mass in the 3 M^ 
model) are excluded by the observations. 

We have measured the effect of rate uncertainties on the two 
most important neutron-producing reactions by independently test- 



ing the high and low rates for '"'C(0!,n)'^0 from [Angulo et al 
(|1999 | NACRE collaborati on) and high, low, and adopted rates for 
^^Ne(ce,n)2^Mg from Uiadis et al. ( 2010 '. Table|4]shows that these 
alternative rates have little effect on sulphur. 

As with most charged particle reactions, the rate of the 
^^Ne(a,n)^^Mg reaction is very highly temperature-sensitive. As 
an extreme example, if the He-burning shell temperature was sm- 
nificantly higher by 33% (0.4 Gl£]instead of 0.3 GK), then the|ll 



|iadis et aTjpOlO ) high rate pred 
3.01 X 10"'' to 1.80 X 10"^ cn 



icts a reaction rate increase from 
cm-'-g^'-s"', i.e., a factor of 631. 
To simulate the effect of a highly increased He-shell temperature 
on neutron production, we test a model with ^^Ne(a,n)^^Mg rates 
boosted by a factor of 600. The surface abundance results of the 
boosted reaction rates are listed in Tableland include: a reduc- 
tion in Ne (0.3 dex), increase in Mg (0.5 dex), an increase in P (0.4 
dex) and a small decrease in S (0.01 dex). The implications of these 
results discussed in Section|4] 



3.3 Comparison to PG1159-035 

The hydrogen-deficient and helium-rich surface chemistry of 
PGl 159 stars is likely caused by a late helium shell flash and conse- 
quent convection zone that extends into the hydrogen-rich surface 
layers l )Schoenbemer|1979[[lberi et al.|1983[|Herwig et al.|1999^ . 
The nucleosynthesis during post-AGB evolution is expected to be 
a relatively small addition to the nucleosynthesis during the pre- 
ceding ~ 20 thermal pulses ( [Herwig et al. 1999), with the excep- 
tion of a very late helium shell flash and proton ingestion episode, 
which may significantly affect the surface composition of light ele- 
ments (up to oxygen) and heavy neutron-capture elements jHerwig| 
|et al.|2011|[Stancliffe et al.|201 \\ . Hence, although our stellar mod- 
els terminate at the tip of the AGB, we expect that our intershell 
abundances of elements heavier than oxygen and much lighter than 
iron will approximately match the surface abundances of PGl 159 
stars. Abundance measurements for elements heavier than oxygen 
in PGl 159 stars are rare in the literature. However, abundances in 
the prototype PGl 159-035 have been measured for many elements 
including Si, P, S, and Fe ( |Jahn et al.|2'007l|Wemer et al.|201 1> and 



where 1 GK = 1 x lO' K. 
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Table 4. Final log e* surface abundances for 3 Mq, Z = 0.0 1 mod el with alternative rate 
sources: '^C(0;,n)'*O low and high estimates from |Angulo et al.|jl999| NACRE co llaboration) 
and ^^Ne(a,n)-'Mg low, median, high, and highx600 rate tables from'Ilia dis et al.|{2o"l Oi. Our 
standard case for comparison uses the ReaclibVO.5 release, which includes fits to the NACRE 
adopted rates for both reactions. 





Ne 


Mg 


S 


P 


CI 


Ar 


Initial 


7.802 


7.472 


6.992 


5.242 


5.062 


6.272 


Standard (JINA ReaclibVO.5) 


8.266 


7.546 


7.024 


5.326 


5.133 


6.302 


'3C(a,n)'^0 NACRE-low 


8.200 


7.537 


7.024 


5.325 


5.135 


6.302 


''C(a,n)'^0 NACRE-high 


8.200 


7.537 


7.024 


5.323 


5.134 


6.303 


22Ne(ce,n)25Mg madis2010-low 


8.207 


7.516 


7.024 


5.319 


5.128 


6.303 


22Ne(ce,n)25Mg madis2010-med 


8.206 


7.517 


7.024 


5.320 


5.130 


6.303 


22Ne(a,n)25Mg Iliadis2010-high 


8.205 


7.517 


7.024 


5.320 


5.130 


6.303 


22Ne(a,n)25Mg Iliadis2010-high-x600 


7.868 


8.068 


7.017 


5.686 


5.146 


6.295 



*\oge(X)=log(Nx/NH) + l2 



Table 5. Top: Observed surface abundances of PGl 159 stars PGl 159-035 and PGl 144+005 with the Sun for comparison. 

Bottom: Intershell abundance results of the models, measured during intershell convection at the last or second last thermal pulse. Rate definitions 
are given in the text. 

References: (l) |Asplund etaLjpOOg) , (2) |Jahn et al.lpOOT) , (3) |Wemer et"aLl(20TT) , (4) |Wemer & Herwig|j2006) 
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Mass Fractions 
















Metallicity Z 


Mass 


C 


O 


F 


Ne 


Si 


P 


S 


Fe 






[M,,] 






[10-7] 


[10--] 


[lO-t] 


[10-5] 


[lO-t] 


[10-4] 


SolarC) 


0.014 


1.00 


0.003 


0.006 


3.66 


0.14 


7.3 


0.58 


3.4 


14 


Scaled Solar 


0.010 




0.002 


0.004 


2.44 


0.09 


4.9 


0.39 


2.3 


9.5 


PG1144+005<'*' 




0.60 


0.570 


0.016 


100 


2.00 










PG1159-035<2-^' 


-Z© 




0.480 


0.170 


32.0 


2.00 


3.6 


0.64 


0.05 


13 



IModels (Z = 0.01) He-Intershell Mass Fractions 



Miniiial 


Rates 


MpMZ 


Mf 


C 


O 


F 


Ne 


Si 


P 


S 


Fe 






[M©] 


[M,.] 






[10-7] 


[10--] 


[lO-t] 


[10-5] 


[lO-t] 


[10-4] 


1.8 


ka02 


2 X 10-3 


0.59 


0.167 


0.012 


202 


2.19 


5.0 


2.2 


2.2 


8.3 


3.0 


kd02 


1 X 10-3 


0.68 


0.176 


0.004 


667 


3.39 


5.2 


1.1 


2.1 


8.8 


3.0 


ka02 


1 X 10-3 


0.68 


0.175 


0.004 


657 


3.39 


5.1 


1.2 


2.2 


8.9 


3.0 


ths8 


NoPMZ 


0.68 


0.183 


0.002 


547 


2.56 


5.3 


5.2 


2.0 


9.4 


3.0 


ths8 


1 X 10-3 


0.68 


0.175 


0.004 


657 


3.39 


5.4 


1.0 


1.9 


9.0 


3.0 


ths8 


5 X 10-3 


0.68 


0.146 


0.006 


809 


6.26 


5.6 


2.4 


1.8 


7.8 


3.0 


ths8 


10 X 10-3 


0.68 


0.123 


0.008 


847 


8.32 


5.9 


2.3 


1.9 


8.0 


3.0 


22Ne-I110-high-x600 


1 X 10-3 


0.68 


0.179 


0.004 


2700 


0.18 


12 


13 


1.6 


1.6 


6.0 


ka02 


No PMZ 


0.98 


0.201 


0.004 


65.1 


1.22 


9.8 


3.5 


1.7 


4.4 



we include them in Table [5] Also shown in Table |5] are the inter- 
shell abundances of our 1.8 Ma,, 3.0 Mrr^, and 6.0 Mr., models for 
comparison. 

I Werner et aT] jl991| estimated the mass of PGl 159-035 at 
0.605^ Q4 Mq by comparison with the evolutionary tracks of 
[Scho enberner (1979). The mass measured by Werner et al.|fl991^ 
is a close match to our 1.8 Mq model, which has a mass at the end 
of the AGB of 0.59 Mq. However, based on more recent evolu- 
tionary calculations by |Miller Bertolami & AlthausH2006^ , |Werner| 
|et al . ( 20 1 1 1 recently revised the mass of PG 1 1 59-059 downwards 
to 0.536+ [J[f) Mq. The predicted final mass of our 1.8 Mq model is 
still contained within the uncertainties of the revised measurement, 
although a better fit may be achieved by a model with a lower initial 
mass. 

In Figure|9] we plot the abundances of PGl 159-035 together 
with the intershell abundances of our 1.8 Mq model. The closely 
matching Fe abundances indicate that our model has a similar ini- 
tial metallicity to PGl 159-035, as low mass stellar evolution leaves 
Fe abundances almost unchanged. We find that our model overpro- 



duces fluorine, with a resulting intershell mass fraction that is 2.0 
times that of PGl 159-035. This may be caused by our model hav- 
ing a higher initial mass than PGl 159-035. 

Our models do not include convective overshoot, and we 
therefore find carbon and oxygen abundances that are too low to 
match the surface of PGl 159-035 ( [Werner et ar i991). The high 
carbon and oxygen abundances of the PGl 159 stars have been 
shown to be consistent with overshoot into the CO core jHerwig| 
|et al.|p^99l l, whereas our abundance results match other standard 
models of AGB stars without convective overshoot jBoothroyd &| 
"Sackmann 1988 ; Kar akas et al.|2002t . Overshoot into the CO core 
could take place during a late helium shell flash, which is when the 
post-AGB star becomes helium-rich, or alternatively might occur 
during every thermal pulse on the AGB. 

We find general agreement between the neon and silicon abun- 
dances, which match to within 10% and 40%, respectively. How- 
ever, our model displays a large overabundance of phosphorus by 
240%, and an even larger overabundance of sulphur, for which our 
model prediction is 44 times too high to match the observation. Our 
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Figure 9. Comparison of our 1.8 M0 intershell abundances with surface 
abundances of PGl 159-035 from |Jahn et aTH2007) and|mrneretal]j20n). 



final intershell sulphur abundances are consistent with the models 
of Herwig that have no extra mixing into the C-O core during ther- 
mal pulses, which predict an abundance of 0.9 times solar jWerner| 
[&Herwig2006j . 



4 DISCUSSION & CONCLUSIONS 

The main finding of this study is that variation in the uncertainties 
that affect nucleosynthesis of AGB stars has little impact on the 
abundance of sulphur in AGB models. The uncertainties that we 
have explored include those associated with the nuclear network 
(e.g., neutron capture cross sections) and those that deal with the 
formation of a '■'C-rich region in the He-intershell of AGB models. 

In terms of the nuclear network, we have measured the effect 
of rate uncertainties of the two most important neutron-producing 
reactions by inde pendently testing th e high and low rates for 
'^C(a,n)'<'0 from Angulo et al.|l |l999| NACRE collaborati on) and 



Iliadis 



high, low, and recommended rates for ^^Ne(a,n)^^Mg from 
|et al.| j20T0j. Table |4] shows that variations in these rates within 
the quoted uncertainties have little effect on sulphur. The rate of 
the ^^Ne(0!,n)^^Mg reaction is highly uncertain, especially at the 
temperatures found in the He-shells of AGB stars (e.g., see discus- 
sions in |Longland et al.|2009]|2012[[Mescher et al.|2012| l. At the 
temperatures of T < 0.30 GK found in the He-shell of low-mass 
AGB stars, this reaction is only marginally activated. However the 
uncertainties quoted at this temperature are considerable, for ex- 
ample, the rate given by Iliadis et al. (2010) varies by a factor of 
1.24 at 0.300 GK, although in comparison the NACRE rate varies 
by a factor of 47 between the upper limit and the recommended 
values at the same temperature. In contrast, the NACRE rate for 
the '■'C(0!,n)'^0 reaction has smaller quoted uncertainties of only 
about 30% at temperatures up to 0.15 GK ( | Angulo et"ar]|1999| 
NACRE collaboration). Note that this reaction is ignited in the He- 
intershell between pulses, when the temperature is much lower than 
during convective thermal pulses. The rate for '■'C(o;,n)'^0 was 
recently re-determined by [Pellegriti et al. ( 2008) a nd Guo et al.| 
(j20T2j. The uncertainties quoted in |Guo et aL j2012^ are now even 
lower, at < 20% at I.O x 10^ K. 

The temperature at the base of the He-flash convection zone is 
dependent on the numerical details of the AGB model and in par- 



ticular on the treatment of the flash-driven convective boundaries 
(e.g., |Herwig|2000l l. The inclusion of convective overshoot at the 
inner boundary of the pulse-driven convective zone (PDCZ) will 
transport additional "^He downward to higher temperatures where 
it burns via the triple- a reaction, leading to increased temperatures 
in the He shefl. In a 3 M©, Z = 0.02 model, |Herwiglj20gO) found 
that the inclusion of diffusive convective overshooting (with over- 
shoot parameter / = 0.016) increased the maximum temperature at 
the base of the PDCZ by 13% from 0.24 GK to 0.27 GK (see also 
Fig. 3 in |Lugaro et al.|2003] l. The application of diffusive overshoot 
to the base of the PDCZ can cause convective mixing to reach into 
the degenerate C-O core, however this phenomenon is not seen in 
3D hydrodynamic simulations of the core-intershell boundary by 
|Stancliffe et al.|j201 1^ . Note that Herwig found modest sulphur de- 
pletion in the He-shell of his model, where the sulphur abundance 
was between 0.6-0.9 times the solar abundance. The sulphur deple- 
tion is directly related to the amount of convective overshoot (and 
therefore presumably the He-shell temperature). For example, Her- 
wig's most recent models calculated with the MESA code l |Paxton| 
let al.|20Il^ and the Nu Grid Multi -zone Post-Processing Network 
tool (Her wig et al.|2008[[Bermett et al.|2012^ and with a lower over- 
shoot parameter applied to the PDCZ of / = 0.008, still show a very 
small sulphur depletion at about 0.8 times the solar abundance. 

Our only reproduction of significant sulphur depletions in the 
He-intershell of our 3M0, Z = 0.01 model resulted from signifi- 
cantly increasing the rate of the ^^Ne(a,n)^^Mg reaction. By boost- 
ing the rate of ^^Ne(a,n)^^Mg by a factor of 600, we obtain a sul- 
phur intershell abundance that is 0.47 times the solar abundance 
(but only ~ 0.7 times the initial abundance) as shown in Table |5] 
As mentioned in Section [X2] the corresponding surface abundance 
predictions relative to our standard rate case include: a reduction in 
Ne (0.3 dex), an increase in Mg (0.5 dex), an increase in P (0.4 dex), 
and only a small decrease in S (0.01 dex). Note that the factor of 
600 is approximately equivalent to a 33% increase in the maximum 
He-shell te mperature from 0.30 GK to 0.40 GK, where the upper- 
limit of the 



Iliadis et al. 



^2010 1 rate increases from 3.01 x 10 " 



to 1.80 X 10^*^ cm^-g^'-s^, i.e., a factor of 631. An increase to 
the He-burning shell temperature of 33% is well beyond the 13% 
increase predicted by Herwig' (20001 with the inclusion of diffusive 
convective overshoot. Furthermore, an increase to the reaction rate 
of this magnitude is well outside of experimental uncertainties, in- 
dicating the difficulty in depleting sulphur via neutron captures in 
AGB models. 

In our investigation of the uncertainties related to '•'C pocket 
formation, we found that the insertion of a partially mixed zone 
(lO^"' M0 and larger) caused a very small reduction in the final 
intershell sulphur abundance by about 5-10% and that the use of 
larger PMZ masses does not necessarily result in greater sulphur 
depletions (Table [Sjl. This is because neutron captures at the top 
of the intershell produce sulphur via phosphorus decay at approxi- 
mately the same rate as they deplete it via decay into chlorine. How- 
ever, because the neutron capture rate of -"'Si decreases as a func- 
tion of temperature while those of phosphorus and sulphur isotopes 
increase, activation of the ^^Ne(0!,n)^^Mg reaction during convec- 
tive thermal pulses at high temperature causes some depletion of 
phosphorus and sulphur. 

Models in the narrow mass range of 2.5 to 3.5 M0 
are known to efficiently produce neon and increase their 
Ne/0 ratio during AGB evolution ^Karakas & Lattanzio||2003] ). 
The elemental increase in neon is due to ^'^Ne production 
from primary '^N in convective pulses via the reaction chain 
''*N(a,7)"*F(/3+)"*0(a,7)22Ne. We find that the inclusion of apar- 
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tial mixing zone significantly increases the final surface abundance 
of neon, and that larger PMZ sizes generally result in higher neon 
abundances (Table [3j. The extent in mass of the PMZ correlates 
with the size of the resulting '^N pocket, which contributes to the 
quantity of '^N available for the production of ^^Ne. 

The positions of our models with large PMZ sizes in the Ar 
versus Ne plane are far from the observational PNe data (Figure 
[Sf and from this we conclude that the partial mixing zone masses 
larger than 5 x lO^-' Mq (50% of the He-intershell mass in the 3 
Mq model) are excluded by the observations. However, this result 
is dependent on the uncertain physics associated with the forma- 
tion of partially mixed zones and '■'C pockets. For example, when 
we inserted the protons into the post-processing code, we made the 
choice that the proton abundance decreases exponentially, i.e., lin- 
early in a logarithmic scale. Studies of the formation of the ^^C 
pocket have found profiles that can slightly differ from this basic 
assumption, as well as from each other (e.g., see discussion in |Lu-| 
Igaro et al.|20T2t . 

To explain the a bundances of young open cluster AGB stars. 



Maiorca et al. 



j2012| find that the effective '"'C required for the ; 



process in low-mass AGB stars is four times larger in models with 
M < 1 .5M0 than that required in more massive AGB stars (but see 
the discussion in [D'Orazi et al.|[20T2l l. |Kamath et al.| ( |2012[ l find 
that PMZ masses of 1.2 x lO^'^ M© or greater improve the fitting 
of AGB models to the C/O ratios and [F/Fe] abundances of stars 
in the Magellanic cluster NGC 1846. Kamath et al.'s best fitting 
PMZ mass represented a fraction of 80% of their 1.86 Mq model's 
1.5 X 10^^ Mp, intershell, with the large PMZ required to produce 
enough fluorine to match the observations. Our results for neon 
and argon in planetary nebulae data are not consistent with such 
a large partial mixing zone. The conflicting results of PMZ stud- 
ies demonstrate the need for further work to identify the formation 
mechanism of '^C pockets in AGB stars. 

In general, AGB models are a good match to fluorine abun- 
dances in post-AGB stars (e.g., [Werner et al.|2005[ l, so the overpro- 
duction of fluorine in our 1.8 M© model relative to PG1159-035 
may be indication that a lower initial mass is needed to model the 
star more accurately. 

Although it was not the original aim of this investigation, our 
finding that neon abundances in PNe can be used to constrain the 
size of the '■'C pocket in low-mass stars is an important and unex- 
pected result. Other investigations of ^^C pocket size in the litera- 
ture exist (e.g., |Bonacic Marinovic et al.|2007| , but our method is 
independent from the uncertainties of AGB s-process abundances. 
To date, the extent and profile of the '■'C pocket are still highly 
uncertain, so the constraint set by neon abundances represents an 
additional clue to understanding this critical step in heavy element 
production via the s-process. A wider study with a finer grid of 
PMZ sizes and more accurate PNe data would likely find a smaller 
upper bound for the PMZ mass and could result in other, new con- 
straints on mixing processes in AGB stars. 

Our results show that variations within the known uncertain- 
ties of nuclear reaction rates and partial mixing zone masses are in- 
sufficient to reproduce the sulphur anomalies in PNe and PGl 159- 
035 via low-mass stellar nucleosynthesis models. We conclude 
that our present knowledge of AGB stellar evolution and the rel- 
evant reaction rates does not support an explanation for the sulphur 
anomaly in terms of the nucleosynthesis in PN-progenitor AGB 
stars. 
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